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Abstract 

We discuss a consistency test which makes it possible to discriminate unknown nu- 
clear background lines from neutrinoless double beta decay with only one isotope. 
By considering both the transition to the ground state and to the first excited + 
state, a sufficiently large detector can reveal if neutrinoless double beta decay or 
some other nuclear physics process is at work. Such a detector could therefore 
simultaneously provide a consistency test for a certain range of Majorana masses 
and be sensitive to lower values of the effective Majorana mass {m v ). 
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1 Introduction 



The first solid evidence for New Physics beyond the Standard Model is the proof that 
neutrinos have a nonvanishing rest mass, as could be convincingly shown in various 
neutrino oscillation experiments [1, 2]. The results show the violation of individual 
flavor lepton number, while conserving total lepton number. A violation of the latter 
would have even deeper consequences for our understanding of the Universe. The gold- 
plated process demonstrating total lepton number violation is neutrinoless double beta 
decay of atomic nuclei. 

Double beta decay is characterized by a change in atomic number Z by two units, 
while leaving the mass number A constant. It is observable for 35 even-even nuclei, as 
single beta decay is energetically forbidden or at least strongly suppressed. It may occur 
in the two-neutrino mode as well as in the neutrinoless mode: 

(Z, A)^(Z + 2, A) + 2e" + 2u e (2u(3(3 decay) (1) 

and 

(Z, A) (Z + 2, A) + 2e~ (Oup/3 decay) . (2) 

2z//3/3 decay conserves lepton number and is therefore allowed in the Standard Model. 
However, this transition is second order in the weak Hamiltonian and thus strongly 
suppressed. Nevertheless, it has been detected experimentally for about a dozen isotopes, 
and their half-lives have been measured. Depending on the nuclide under consideration, 
half-lives are of the order 10 20 y. 

The second process, Ozy/5/3 decay, is forbidden in the Standard Model as it violates 
total lepton number by two units. It can be seen as two subsequent steps ("Racah 
sequence" [3]): 

(Z,A)^(Z+l,A) + e- + v e 
{Z + 1, A) + v e (Z + 2, A) + e~ , U 

which is only possible if the neutrino is its own antiparticle (a so-called Majorana par- 
ticle) and if it has a nonvanishing rest mass to account for the helicity matching. As 
the decay is a kind of a black box, any AL = 2 process could contribute to 0v(3(3 decay. 
OvfiP decay has not been observed experimentally yet. 1 Nevertheless, stringent bounds 
on the half-lives of various elements have been extracted from experiments. Best limits 
are of the order 10 25 y. 

The total decay rate of neutrinoless double beta decay (if mediated by light neutrino 
exchange) is given by [6] 



m u f M^'G [,U (Q,Z). (4) 



in 2 r°; 2 



2 



1 Note that one claim for a positive signal of 0^/3/3 decay exists. A subgroup of the Heidelberg- 
Moscow collaboration gives the half-life T®" 2 = 1.98 x 10 25 y, where the 3c range is given by (1.04- 

20.38) x 10 25 y [4], and a more recent analysis results in T°Y 2 = (2.23±g;|i) x 10 25 y [5]. 
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Figure 1: Schematic plot of the sum energy spectrum of the emitted electrons for the 
various decays to be discussed in this paper. As emitted neutrinos may take away an 
arbitrary amount of energy, the spectrum for 2u(3/3 decay is continuous, whereas the 
spectrum of 0v(3(3 decay is a single line at the maximum energy Q. The ratio between 
both is not to scale. For the decay to an excited final state, the diagrams are qualitatively 
the same. However, as energy is taken away by the emitted photons, the line for 
decay lies at lower energies. Moreover, the number of decays to excited states is lower 
than the number of decays to the ground state. 



Here, T®Y 2 is the half-life of Ov/3/3. The nuclear matrix element M° v and the phase 
space integral G 0l/ (Q, Z) depend on the nucleus under consideration. In 0v(3(3 decay, it 
is thus possible to measure (m Ue ), the so-called effective Majorana mass of the electron 
neutrino. For light neutrinos, it is given by 



m, 



1 2 e 2ia -m, 



(5) 



where U e i are the elements of the first line of the Pontecorvo-Maki-Nakagawa-Sakata 
neutrino mixing matrix, and a,i are the two Majorana CP-violating phases. In case of 
CP-invariance (ctj = 0, n), we may write 



mxUly ± m 2 Ul 2 ± m 3 U^ 3 



(6) 



The experimental signal of 0v(3(3 decay is two electrons in the final state, whose 
energies add up to the Q value of the nuclear transition. For 2v(3(3 decay, however, 
the sum energy spectrum of both electrons will be continuous as the neutrinos may 
take away an arbitrary amount of energy (see Fig. 1). The total decay rates, and 
hence the inverse half-lives, are a strong function of the available Q value. The rate 
of Qvf3f3 decay scales with Q 5 compared to a Q u dependence for 2v(3f3 decay, due to 
the dependence on the phase space factor G 0l/ (Q,Z) [cf. Eq. (4)]. Therefore, isotopes 
with a high Q value (above about 2 MeV) are usually considered for experiments on 
0v(3(3 decay. This restricts the candidates to 11 promising isotopes, which are given in 
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Figure 2: Scheme of the double beta decay of 76 Ge. Single beta decay to 76 As is en- 
ergetically forbidden, thus the only open decay channel is double beta decay. It may 
proceed in the two-neutrino mode, as well as in the neutrinoless mode. Additionally, we 
will discuss the decay to the first excited + level. 



Table 1 together with the corresponding Q values. For recent reviews on double beta 
decay, see [6-9]. 

2 Double beta decay to excited 0^ states 

Usually, double beta decay to the ground state ("g.s." in equations) of the final nucleus 
is considered. However, practically all interesting nuclei, i.e., those with a Q value above 
2 MeV (see Table 1) have at least one excited + and one excited 2 + state which is 
accessible by double beta decay as well. The level scheme of 76 Ge is given in Fig. 2 as an 
example. Transitions to excited 2 + states might be dominated by potential contributions 
of V+A interactions (see, however, [20]). 

The decay rate to excited states is lower due to the lower Q value of the decay. The 
ratio between the decay rate to the excited Of state and the ground state is given by 

ry _ (Q-E(of)r (M°t\ 2 

r g . s . Q n \m* b -) ' 1 ' 

where E{0f ) is the energy of the first excited state with respect to the ground state. 

For 0v(3(3 decay, n = 5 describes the phase space dependence, A4q^' denotes the nuclear 

+ 

matrix element for the decay to the ground state, whereas M.q v denotes the nuclear 
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Table 1: All double beta emitters with a Q value larger than 2 MeV and their corresponding Q values. All Q values with an 
error larger than 1 keV are taken from [10]; all others were recently remeasured using Penning traps [11-15]. Furthermore, 
the energy of the first excited + state as taken from [16] is shown. All nuclear matrix elements are obtained within the 
IBM-2 model [17-19]. As can be seen from pure phase space considerations, 150 Nd and 96 Zr would be the choice; however, 
including the matrix elements, 150 Nd and 76 Ge seem to be the most promising cases to study. 



Decay mode 


Q [keV] 


E(0f) [keV] 






(Q - £(0+)) 5 /Q 5 




r 0+/ r g-s- 




zu zz 


4274 ± 4 [10] 

L J 


2997 






2.38 x 


10~ 3 








76 Ge _v76g 


203Q 04 ± 16 [111 


1122 


5 465 


2 479 


1.84 x 


io- 2 


206 

\J . Zj \J\J 


3 79 x 10" 


-3 


82q p v82 K 
34 oe— i-ggivi 


2995.5 ± 1.9 [10] 


1488 


4.412 


1.247 


3.23 x 


io- 2 


0.080 


2.58 x 10' 


-3 


$Zr-^Mo 


3347.7 ±2.2 [10] 


1148 


2.530 


0.044 


1.22 x 


IO' 1 


3.02 x 10" 4 


3.70 x 10' 


-5 


^Mo-^Ru 


3034.40 ±0.17 [11] 


1130 


3.732 


0.419 


9.74 x 


io- 2 


1.26 x 10~ 2 


1.23 x 10" 


-3 


JS-Pd^Cd 


2004 ± 11 [10] 


1473 


3.623 


1.599 


1.31 x 


IO" 3 


0.195 


2.54 x 10" 


-4 


^Cd-^Sn 


2809 ± 4 [10] 


1757 


2.782 


1.047 


7.37 x 


IO" 3 


0.142 


1.04 x 10" 


-3 


^Sn^Te 


2287.8 ± 1.5 [10] 


1657 


3.532 


2.721 


1.59 x 


10~ 3 


0.594 


9.46 x 10' 


-1 




2527.518 ±0.013 [13] 


1794 


4.059 


3.090 


2.06 x 


10~ 3 


0.580 


1.19 x 10' 


-3 




2457.83 ± 0.37 [14, 15] 


1579 


3.352 


1.837 


5.84 x 


10~ 3 


0.300 


1.76 x 10" 


-3 


^Nd-^Sm 


3371.38 ±0.20 [12] 


740 


2.321 


0.395 


2.90 x 


IO' 1 


2.90 x IO" 2 


8.39 x 10- 


-3 



matrix element for the decay to the first excited + state. A similar relation holds for 

the 2is/3/3 decay mode, however with different matrix elements M. 2 l> an d -M-fu as wen 
as with a different scaling law for the phase space, i.e., n — 11. 

So far, excited state transitions have only been observed in two nuclides, namely 
100 Mo [21] and 150 Nd [22], both considered to be 2i>/3f3 decay. The question arises 
whether, in case of an observation of the ground state transition within a single iso- 
tope, the transition to the first excited + state can be used as a consistency check, 
i.e., whether it is possible to prove 0^/3/3 decay with one isotope in one experiment in 
two different ways. Such a test might be desirable in future large-scale experiments 
due to the involved costs. The results for all double beta emitters with a Q value of 
at least 2 MeV are compiled in Table 1. For consistency, the matrix elements from a 
single method, the interacting boson model (IBM-2), have been used [17-19]. Thus, 
numbers might change slightly if other calculations are used. Unfortunately, there is no 
complete set of such matrix elements available for all 11 isotopes, including also excited 
state transitions. It turns out that the two most suitable choices for such an internal 
consistency check would be 150 Nd and 76 Ge, the first one about a factor of 2.21 better. 

Let us remark that a possible benefit of the proposed method is that the nuclear 
matrix elements of the transition to the ground state and to the excited state may have 
common uncertainties, which would cancel in the ratio. Improvement of the nuclear 
matrix elements will allow for a more precise extraction of the effective neutrino mass 
from half-life measurements. 



3 Experimental considerations 

The expected signature for the required second decay mode into the excited Of state 
will be two electrons and two gammas with defined energies in contrast to the ground 
state transition having only two electrons with a defined sum energy. Thus, the gammas 
must be clearly separated from the emitted electrons in the experiment, otherwise they 
would look like a ground state transition. Hence, a purely calorimetric approach without 
spatial resolution to determine the individual gammas will fail. Consequently, in a large 
homogeneous detector, there must be spatial resolution to see the gammas independently 
from the emitted electrons. In a high granularity detector, the granularity should be 
chosen such that both gammas might leave the crystal containing the decay without any 
interaction, making it possible to search for coincidences with high efficiency. 

As all double beta decays into the first excited + state will de-excite via the sequence 
of + — > 2 + — > + , there will also be a 77 angular correlation, which for the given angular 
momentum sequence is 

W{9) = | x (1 -3cos 2 # + 4cos 4 #). (8) 
8 

This function is plotted in Fig. 3. It is easy to see that the angles and ir have the 
highest probability. 
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Figure 3: Angular correlation from Eq. (8). 

Which types of background are to be expected? In the chart of nuclides, double 
beta emitters are surrounded by unstable isotopes. Thus, the intermediate nucleus in 
the double beta system of interest — which also might be produced by (p, n) reactions on 
the double beta emitter — is unstable and its beta decay into excited states will lead to 
the same gamma- signature. However, the energy spectrum of the single beta decay will 
be continuous but overlapping with the double beta electron signal. The fraction of the 
beta electrons in the peak range depends on the energy difference of the ground states of 
the double beta emitter and the intermediate nucleus. If it is small, only electrons close 
to the endpoint of the beta spectrum will contribute. If it is large, more electrons will 
contribute. A detailed estimate depends also on the quantum numbers of the ground 
state of the intermediate nucleus, as allowed or forbidden beta decays will lead to different 
electron energy spectra. Thus, it is essential to measure the electron energy accurately 
or to build a detector which is able to discriminate one and two electrons, typically done 
in detectors with tracking capabilities. 

Let us comment on the possible origin of the aforementioned (p, n) reactions. Nor- 
mally, there are no free protons in an underground experiment, so these reactions are not 
an issue. But as our consistency test relies on a very low background, even tiny contri- 
butions have to be considered. Such a contribution is the production of protons by (n,p) 
reactions on a nucleus. Underground, high energy neutrons are dominantly produced 
by muon interactions in or close to the experiment. These neutrons in principle have 
enough energy to do (n,p) reactions. A detailed estimation, however, depends on the 
actual cross section for (n,p) reactions (mb region for neutron energies below 100 MeV) 
and the following (p,n) reactions (mb— b region, depending on the proton energy). 

Also, external backgrounds may be an issue. The signal of a decay into excited states, 
however, will be a triple coincidence with well-defined energies of all involved particles. 
Additionally, angular correlations exist (at least between the gammas), and the total 
sum of particle energies must correspond to the Q value of the double beta decay. These 
constraints make the signal search more or less background free, of course depending a 
little bit on the detector technology used. Especially in the case of Ge detectors with 
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their superb energy resolution, the triple coincidence is so sharply defined that it cannot 
be mimicked by any other process. Decay sequences with the same gamma energies are 
very unlikely, furthermore the electron in the sequence will be continuous, and only a 
small fraction will have the right energy. Also, triple Compton events (note that we have 
three different energy depositions) are very unlikely. Moreover, applying the equation of 
Compton scattering to the three energy depositions will immediately tell whether this 
is consistent with Compton scattering or not. 

The major background will be the 2z//3/3 decay into the Of excited state, observed 
so far for two isotopes. As the high energy tail of the continuous two electron spectrum 
can mimic the signal, energy resolution becomes the crucial experimental quantity for 
rejecting this background. 

To be more specific, the two most promising nuclides are discussed in more detail. 
First, consider the case of 76 Ge which is used in the GERDA and MAJORANA exper- 
iments. In addition to the signal in the sum energy spectrum of the electrons in form 
of a peak at 2039 keV as an indication of the ground state transition, the decay into 
the first excited + state would be indicated by a sum energy of the electrons of 917 
keV associated with two gammas of 559.1 keV and 563.2 keV, respectively. Clearly, a 
coincidence measurement would be preferential for this form of decay. Even ignoring 
the angular correlation among the photons, this channel should be observable more or 
less free of background. In this case, one or two events would indicate an observation, 
which, however, implies 264-528 neutrinoless ground state transitions (see rates in Ta- 
ble 1). First Monte Carlo simulations show that for 76 Ge in the case of using detectors 
in form of disks of 15 cm diameter and 1 cm thickness, about 60 % of the gammas 
are expected to leave the crystal without interaction, future simulations might optimize 
this number. The photons might be detected in neighboring Ge detectors or an active 
medium surrounding the crystals. 

The most promising candidate is 150 Nd, which is considered to be used in SNO+ 
and DCBA. Moreover, it is an option for SuperNEMO. It has a Q value of 3371 keV 
with two gammas of 334 keV and 406 keV, respectively. As before, the independent 
measurement of all these energies requires high granularity detectors, relatively small 
crystals in a liquid with a fair spatial resolution or tracking devices. Again, in case of 
1-2 excited state events, this would imply 120-240 events in the neutrinoless ground 
state transition. In general, depending on the actual value of the effective Majorana 
mass, this can be a challenging measurement. 

Hence, how realistic is such an approach? Can the necessary number of counts 
(several hundred, see above) be reached for the decay into the ground state in future 
detectors? 

In case the half-life T\i<i of the isotope under consideration is much longer than the 
measuring time t, we may write the number of double beta decays as 

%=^, (9) 

il/2 

where a is the isotopical abundance of the nuclide of interest, M is the used mass and 
Na is the Avogadro constant. However, in experiments, we may be confronted with 
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Table 2: Running times to be accumulated for the two types of next-generation detec- 
tors discussed in this paper, to be able to use the proposed consistency test. The two 
thinkable scenarios are described in the text. 







76 Ge 




i50 Nd 


Klapdor's claim Ty 2 = 2.23 x 10 25 y [5]: 
background free 
background limited 


Mt 
Mt 


= 2.0 tony 
= 121 tony 


Mt 
Mt = 


= 3.4 tony 
= 30.9ktony 


(m„J = 50 meV: 
background free 
background limited 


Mt 
Mt = 


= 48.3 tony 
= 69.9ktony 


Mt 
Mt -- 


= 14.2 tony 
= 543 kton y 



background, such that there are two different possible dependencies of the expected 
half-life sensitivity: 

(^1/2) 1 aMet (background free) (10) 

or 



(^1/2 ) 1 ae \ n> I ?i (background limited) . (11) 
V B £\E 

Here, e is the efficiency for detection, B is the background index [typically quoted in 
counts/(keV kg y)], and AE is the energy resolution at the peak position. See [23] for 
a more detailed discussion on 0^/3/3 experiments. 

Consider two types of next generation detectors: a 76 Ge detector, and a 150 Nd de- 
tector. Two scenarios are thinkable: the Klapdor claim TfY 2 = 2.23 x 10 25 y [5] is right, 
so one only would have to reach this half-life. However, it is not improbable that the 
effective Majorana neutrino mass is as low as 50 meV, where the inverted hierarchy be- 
gins (cf. Fig. 4). Results for the running times to be accumulated in both scenarios to 
be able to use the proposed consistency test are given in Table 2. All parameters used 
are given in Table 3. 

How big may the half-life maximally be, so that the proposed consistency check can 
be used in a running time of, say, 10 years? This means that we need at least 264 or 120 
decays to the ground state in 76 Ge or 150 Nd, respectively. Assuming a background-free 
experiment, and using the values from Table 3, for a 76 Ge detector of 1 ton, the maximal 
half-life where this number of decays to the ground state is reached is 

Ti /2 = 1.11 x 10 26 y- (12) 

This corresponds (using the matrix elements and phase space factors in Table 3) to an 
effective Majorana neutrino mass of 0.11 eV. For a 150 Nd detector of 200 kg, the maximal 
half-life is 

T 1/2 = 1.31xl0 25 y, (13) 
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Table 3: Parameters used in the calculations. The values for 76 Ge are as currently 
taken for the GERDA experiment (phase I) [24], the values for 150 Nd are as expected 
for SuperNEMO [25]. Nuclear matrix elements are IBM-2 [17-19], and the phase space 
factors are taken from [26]. 





76 Gp 


150 NH 


t -i- ■ i u j r —li 
Isotopicai abundance a [m u J 


nrQ/ ro A 1 

85% [24J 


90% 


hvmciency e 


0.60 


0.30 [25] 


Number of decays to ground state Npp 


264 


120 


Background B [counts/ (keV kg y)] 


0.01 [24] 


0.02 [25] 


Energy resolution at peak position AE 


3 keV [24] 


4% [25] 


M 0u 


5.465 


2.321 


G 0v [y- 1 eV" 2 ] 


0.25 x 10~ 25 


8.03 x 10~ 25 


typical mass [kg] 


1000 


200 



which corresponds (using again Table 3) to an effective Majorana neutrino mass of 
0.14 eV. Thus, for detectors of the size of several hundred kilograms up to lton, the 
method should work down to about 100 meV. 

4 Summary and conclusions 

We have proposed a method to check within a single experiment whether a possibly 
observed signal in a future Oz//3/3 detector is really due to Oz//3/3 or due to some unknown 
nuclear line. This question will arise if a positive signal is seen in a second-generation 
detector like GERDA, which is about to start data-taking and which will then reach and 
exceed the sensitivity of the Heidelberg- Moscow experiment to a certain point. Usually, 
it is argued that another second generation experiment with a different isotope can settle 
the question. We pointed out that it is also possible to combine effort into one large 
detector, instead of using various different isotopes. Such a detector would therefore 
serve different purposes: It would have sensitivity to lower (m„ e ) and it would be able 
to check a claim for higher (ra„ e ) due to the very characteristic features described in 
this paper. See Fig. 4 for the reach for both modes. It is clear from the figure that 
the proposed consistency test is only viable for relatively large effective neutrino mass. 
Should we really have to cover the whole inverted hierarchy, or even go down to normal 
hierarchy, one large detector may not be viable due to the large amount of material 
needed, and a set of various isotopes in smaller detectors may be preferable to check for 
consistency. We mentioned also that our proposed method may allow a better extraction 
of the effective neutrino mass, since common errors in the calculations of nuclear matrix 
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Figure 4: The effective Majorana mass (m Ue ) as a function of the lightest neutrino mass 
eigenvalue for inverted (blue/light blue) and normal (orange/yellow) hierarchy. In the 
quasi degenerate case, both hierarchies overlap. Bold colors denote the best fit values 
range (varying the CP phases), light colors give the corresponding 3a values range. 
The best fit value (m Ue ) = 0.34 eV obtained in the Heidelberg- Moscow experiment [4] is 
marked. A future 1 ton 76 Ge experiment [27] could probe the inverted hierarchy down 
to (m Ue ) = 0.01 eV. Using this experiment for the consistency test, the probed value for 
(m Ve ) would be higher, due to the lower rate of the decay to excited states. 



elements may cancel in ratios. Of course, a deeper discussion of the consistency test 
is necessary before a particular experimental setup may be proposed. However, our 
considerations should have proven the viability of such a method. 
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